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Quark-nuclear matter (QNM) is a many-body system containing hadrons and 
deconfined quarks. Starting from a microscopic quark-meson coupling (QMC) 
Hamiltonian with a density dependent quark-quark interaction, an effective quark- 
hadron Hamiltonian is constructed via a mapping procedure. The mapping is 
implemented with a unitary operator such that composites are redescribed by 
elementary-particle field operators that satisfy canonical commutation relations in 
an extended Fock space. Application of the unitary operator to the microscopic 
Hamiltonian leads to effective, hermitian operators that have a clear physical inter- 
pretation. At sufficiently high densities, the effective Hamiltonian contains interac- 
tions that lead to quark deconfinement. The equation of state of QNM is obtained 
using standard many-body techniques with the effective quark-hadron Hamilto- 
nian. At low densities, the model is equivalent to a QMC model with confined 
quarks. Beyond a critical density, when quarks start to deconfine, the equation of 
state predicted for QNM is softer than the QMC equation of state with confined 
quarks. 



1. Introduction 

One of the most exciting open questions in the study of high density hadronic 
matter is the identification of the appropriate degrees of freedom to describe 
the different matter phases. For systems with matter densities several orders 
of magnitude larger than the nuclear saturation density, one expects a phase of 
deconfined matter composed of quarks and gluons whose properties very likely 
can be described by perturbative QCD. For ground state nuclei, there is a large 
body of experimental evidence that their gross properties can be described 
more economically employing hadronic degrees of freedom, rather than quarks 
and gluons. On the other hand, for matter at densities not asymptotically 
higher than the saturation density, like the ones in dense stars and produced 
in high-energy nuclear collisions, the situation seems to be very complicated, 
since hadrons and deconfined quarks and gluons can be simultaneously present 
in the system. Presently, it is not possible to employ QCD directly to study 
such systems and the use of effective, tractable models are essential for making 
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progress in the field. 

One attractive model to study the different phases of hadronic matter in 
terms of explicit quark-gluon degrees of freedom is the quark-meson coupling 
(QMC) model, originally proposed by Guichon and subsequently improved by 
Saito and Thomas 1 . For a list of references on further improvements of the 
model and recent work, see Ref. 2 . In the QMC model, matter at low density is 
described as a system of nonoverlapping MIT bags interacting through effective 
scalar- and vector-meson degrees of freedom. The effective mesonic degrees of 
freedom couple directly to the quarks in the interior of the baryons. At very 
high density and/or temperature, when one expects that baryons and mesons 
dissolve, the entire system of quarks and gluons becomes confined within a 
single, big MIT bag. 

In a regime of very high density, the description of hadronic matter in terms 
of nonoverlapping bags should of course break down, since once the relative 
distance between two bags becomes much smaller than the diameter of a bag, 
the individual bags loose their identity. The density for which this starts to 
happen is presently unknown within QCD. 

In the present communication we introduce a generalization of the QMC 
model that allows to include quark deconfincment at high density. Our start- 
ing point is a relativistic quark potential model 3 . From the model quark 
Hamiltonian, we construct a unitarily equivalent Hamiltonian that contains 
quark and hadron degrees of freedom. Starting from the Fock-space represen- 
tation of single-hadron states, a unitary transformation is constructed such that 
the composite-hadron field operators are redescribed in terms of elementary- 
particle field operators in an extended Fock space. When the unitary transfor- 
mation is applied to the quark Hamiltonian, effective, hermitian Hamiltonians 
with a clear physical interpretation are obtained 4 . In particular, one of such 
effective Hamiltonians describes the deconfinement of quarks from the inte- 
rior of the hadrons. The equation of state of QNM can be calculated using 
standard many-body techniques with the quark-hadron Hamiltonian. We will 
show that at low densities, the model is equivalent to the QMC model and, 
beyond a critical density, when quarks start to deconfine, the equation of state 
predicted for QNM is softer than the QMC equation of state. 



2. QMC model with confined quarks 

The nucleons are bound states of three constituent quarks. Constant scalar 
(<7o) and vector (w ) meson fields couple to the constituent quarks in the inte- 
rior of the nucleons. Each constituent quark satisfies a Dirac equation of the 



3 

form 

-iS-f + f3°m* + 1/2(1 + 0)V(r) 

where 

™* =m q -g%a , E* = e* - g% w 0) V(r)=ar. (2) 

The only difference with the model of Toki et al. 3 is the form of the potential, 
while theirs is a harmonic oscillator, ours is a linearly rising one. For a linearly 
rising potential, the Dirac equation cannot be solved analytically. We use 
the saddle point variational principle (SPVP) 5 to obtain an approximated 
solution. Since the Dirac Hamiltonian does not have a lower bound for the 
energy, the traditional nonrelativistic variational method cannot be employed. 
The SPVP amounts to minimizing (maximizing) the energy expectation value 
with respect to the variational parameters corresponding to the upper (lower) 
component of the Dirac wave function. We use as ansatz for the Dirac wave 
function 5 

tf(r)= ( - U{ ~ \ Axs, (3) 
\ia ■ rv(r) ) 

with 

u(r) = Ne~ >?r ' 2 ' 2 , v{r)=ij/Xa-Wu(r), (4) 

where N is a normalization constant, and A and 7 are the variational param- 
eters. The parameters are found by minimizing the energy eigenvalue e with 
respect to A and maximizing it with respect to 7. We note here that for a 
harmonic oscillator potential, the SPVP with this ansatz leads to the exact 
solution. 

Following the traditional path in the QMC model, we initially fix the pa- 
rameters of the model in vacuum. The nucleon mass in vacuum (er = = wo) 
is given by 

Mjv = 3£-£ A, (5) 

where the last term above is used to take into account the cm. energy of 
the three-quark state and other short-distance effects not taken into account 
by the confining potential, such as gluon exchange. Sq is fitted to obtain 
M N = 939 MeV. The value of the string tension is taken to be a = 0.203 GeV 2 
and m q = 313 MeV. With these parameters, the SPVP leads to A = 2.38 fm _1 
and 7 = 0.346. The value required for eo to fit the nucleon mass is 4.67 MeV fm. 

Next, we proceed to obtain the energy of nuclear matter. Nuclear matter 
in the QMC model is modeled as a system of nucleons treated in the mean 
field approximation, in which the quarks remain confined within the nucleons. 



V>(0 = E* q ^{r), (1) 
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The nucleon mass is now obtained as above, but now including the mean fields 
coupled to the quarks. The energy density of symmetrical nuclear matter is 
given by the traditional expression in the QMC model 

E f kp d 3 v 1 1 

V = 4 Jo E * N ^ +39 ^ PB + 2 m2a<j2 ° ~ 2 m ' W °' (6) 

where = p 2 + Mff, ps is the baryon density and m a and m w are the 
masses of the mesonic excitations. The next step consists in determining the 
mean fields. The vector mean field, from its equation of motion, is simply given 
in terms of pb, and the scalar field is obtained by minimizing E with respect 
to <to, as usual. The coupling constants are then obtained by fitting E to the 
binding energy of nuclear matter at the saturation density, i.e. E/B — Mjy = 
-15.7 McV at p B = po = 0.17 fin" 3 (or k F = 1.36 fin" 1 ), where B is the 
baryon number - in this case B is equal do the number of nucleons. Note that 
for each value of pb, one has to use the SPVP to obtain the in- medium values 
of A and 7. The values obtained for the coupling constants are g% — 6.1355 
and ig^j — 6.285. The incompressibility is found to be K — 248 MeV. 

In the next section we generalize the model to allow the deconfincment of 
quarks. 

3. QMC model with quark deconfinement 

Here we construct an effective Hamiltonian that contains hadron and quark 
degrees of freedom. The starting point is the quark model discussed in the 
previous section. In this model, the one-nucleon state can be written in a 
second-quantized notation as 

\a) = fit |0), Bl = ^£™ 9 t i3 t 29 t 3) (7) 

where the q^s are constituent-quark creation operators and -fy^ 1 ^ 2 ^ [ s ^ e 
Fock-space nucleon amplitude - for independent quarks, this is simply the 
product of three single-quark wave functions. The convention of summing 
over repeated indices is used throughly. The quark creation and annihilation 
operators satisfy the usual canonical anticommutation relations 

{%,<&} = 5^, {q»,qv} = {ql,qt} = o. (8) 

The index a denotes the spatial and internal quantum numbers, such as inter- 
nal and cm. energies and the spin-isospin quantum numbers of the nucleon. 
Similarly, the quark indices p, identify the spatial and internal quantum num- 
bers as momentum, spin, flavor and color. The amplitude vf/w^Ms j s taken to 
be orthonormalized: 

(a\(3) = $*fiWM3^MiWM = Saf3 . (9) 
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In the abbreviated notation we are using, the Hamiltonian corresponding 
to Eq. (1) can be written as 

H q =T q + V qq = T (n) qfa + ^V qq {jiv, op) q\qlq p q„, (10) 

where V qq is the confining potential. 

Using the quark anticommutation relations of Eq. (8) and the normalization 
condition of Eq. (9), one can shown that the nucleon operators, B a and B^, 
satisfy the following anticommutation relations 

{B a ,Bl} = 5 aP -A aP , {B a ,B } = O, (11) 

where 

A a p = M* a m3 % m <il 3 % 3 le^^J'^iifefc. (12) 

In addition, one has 

{%,bD = yf *™<i<4, {Q^ B «} = °- ( 13 ) 

The term A a p is responsible for the noncanonical nature of the baryon anti- 
commutator. This term and the nonzero value of Eq. (13) are manifestations 
of the composite nature of the baryons and the kinematical dependence of the 
quark operator and nucleon operators. This fact complicates enormously the 
mathematical treatment of many-body systems in which deconfined quarks 
and nucleons are simultaneously present. The mapping formalism of Ref. , 
known as the Fock-Tani (FT) representation 6 , is a way to circumvent such 
complications. We will shortly review this formalism in the context of the 
present model. For further details, and applications for more general models, 
the reader is referred to Ref. 4 . 

The basic idea is to extend the original Fock space by introducing fictitious, 
or ideal nucleons that satisfy canonical anticommutation relations. The unitary 
operator is constructed in the extended Fock space such that 

|a> = st|o>^t/->) = |a) = ^|o), (u) 

where ideal baryon operators w a and b a satisfy, by definition, canonical anti- 
commutation relations 

= <W> {&«,&/»} = 0. (15) 

The state |0) is the vacuum of both q and b degrees of freedom in the new 
representation. In addition, in the new representation, the quark operators q^ 
and q are kinematically independent of the b^ a and b a 

{qM = {q^,bl t } = 0. (16) 
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The unitary operator U can be constructed as a power series in the bound 
state amplitude The rational for this is clear: in situations that the quarks 
remain confined in the interior of the nucleons, the term A Q( 3 plays no role, can 
be taken to be zero and the unitary operator becomes trivial 6 4 . This is the 
situation for low densities, when the internal structures of the nucleons do not 
overlap significantly in the system. As the density of the system increases, the 
quark structures of different nucleons start to overlap. An expansion in powers 
of "f's offers a power counting procedure to construct the unitary operator. 

The effective Hamiltonian is constructed by applying the unitary operator 
to the microscopic quark Hamiltonian of Eq. (10), H e ff = U~ 1 H q U. The 
zeroth-order U is trivial and not interesting. The first-order U brings interest- 
ing effects. At this order, we denote the effective Hamiltonian by H^j, where 
the superscript (1) means that U has been evaluated up to the first order in 
^. H^jf can be written as 

H% = T q + H b + V qq + V qb + ■ ■ ■ . (17) 

The • • • refer to terms not relevant for our discussion here. Hb = Tb + Vbb, where 
Tf> is a single-nucleon energy and Vbb is an effective nucleon-nucleon interaction 
without quark exchange. This term leads to the normal QMC model, in which 
the many-body system is described by nonoverlapping nucleons - no quark- 
exchange. In particular, it can describe Fock terms in the QMC model 7 . The 
term V qq contains two-quark and three-quark interactions. It can be shown 
that if ^ is a bound-state eigenstate of the original quark Hamiltonian, V qq 
collapses to 

V qq = ^ V qq (fiv; op) q\q\q p q a - E a B\B a . (18) 

It is not difficult to show that this interaction leads to a quark Hamiltonian 
that has a positive semidefinite spectrum. That is, after the transformation, 
the resulting Hamiltonian involving only quark operators is unable to bind 
three quarks to form a nucleon, it describes only states in the continuum. The 
term V qb is given by 

I 

- H(^;ap)^ pT3 A(^f, 2 ^; ^r 3 )l ^^g^ + h.c. , (19) 



Vbq = — 



if(/ii/z 2 ;CTp)* 



O-Pfj.3 



where h.c. denotes hermitian conjugation and A(pvT;ap\) — J2 a ^^ T ^*^ pX 
is known as the bound-state kernel. If * is a stationary state of the microscopic 
quark Hamiltonian, one obtains 

V bq = 0, (20) 



7 



since A(/iz/t; apX) 1 ^^ — . This result reflects the stability of the baryon 
bound state to spontaneous decay in the absence of external perturbations. 
This is clearly the case for a nucleon in vacuum. Also, in the QMC model, when 
one constructs the unitary transformation U with a IP that is an eigenstate of 
the microscopic quark Hamiltonian with the mean fields Co and u>o m it) the 
term Vb q continues to be zero, and explicit quark degrees of freedom will not 
be present in the system at this order of ^ . 

Now, in a many-body system, the confining quark-quark interaction will 
become modified due to a variety of effects. Some of such effects, as self- 
energy corrections from quark loops, can be calculated within the model using 
standard many-body techniques. However, in a high density system there are 
other QCD effects that are not captured by the model, such as pair creation 
and gluonic interactions, that eventually will lead to quark deconfinement. 

The formalism we just described allows to include in an effective way de- 
confinement in the QMC model. One generates an effective quark-hadron 
Hamiltonian as above using 'J's that are eigenstates of the QMC Hamiltonian, 
Eq. (1). Now, if V(r) is modified such as that it does not lead to absolute 
confinement, the term Vb q , given by Eq. (19), is not zero. In a mean field per- 
spective, the Hamiltonian of Eq. (17) leads to two Fermi seas, one for baryons 
and one for quarks. The crucial, and difficult point here is to obtain the rela- 
tive abundances of baryons and quarks in the system. This can be evaluated 
in an approximated way as follows. 

Let Z be the fraction of baryons in the system, 

Y J (blb a )=ZB, 5>tg„> = (1 - Z)B, (21) 

where B is as in the previous section the total baryon number. In the mean field 
approximation - or independent-particle approximation - and for sufficiently 
small Vbq, Z can be estimated by the perturbative formula 

Z-' = l + {b\V^ ql ^-^Vb q \b), (22) 

where H a is T q + T b , and P = 1 — \b)(b\ is a projection operator. 

In order to evaluate Eq. (22), we postulate a density dependence for the 
confining interaction of the form 8 

V(r)=ore-^ r \ (23) 

where /j, is a prescribed function of ps- We use a simple formula for ji, such 
that it is zero for baryon densities below three times the normal nuclear matter 
density p 0} and for higher densities it increases linearly with the density as 
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r(fm) 

Figure 1. The confining potential in vacuum (solid line) and in matter for two different 
baryon densities. 

p, = pb/'^Po — 1- In Fig- (1) we show the potential of Eq. (23) for zero density, 
and 5 and 10 times the saturation density of normal nuclear matter. 
The energy density of the system can be written as 

E f F d 3 p 1 2 2 1 2 2 

V = J (2^)3 n{ - P > + 9ujUJ ° PB + 2 mffCT ° ~ 2 mujUJ ° 

+ 12 1 j^E^k), (24) 

where E*(k) = ^Jk 2 + to* 2 and the Fermi momenta k b F and k q F are related to 
the nucleon density and quark density as 

Pb = Z PB = 2(4) 2 /3tt 2 , p q = (1 - Z)pb = 2(4) 2 /tt 2 , (25) 

At this point, it is important to notice that since p only starts to operate for 
densities larger than three times the normal density, the coupling constants 

and g% are the same as before. Of course, for higher densities, there is a 
somewhat complicated self-consistency problem to be solved, since Z is density 
dependent. Therefore, in the process of obtaining a, the iterative problem 
becomes more complicated. 

In order to proceed, we need Z as a function of ps- It can be calculated 
numerically with our ansatz wave function given above. The calculation, how- 
ever, involves multidimensional integrals that must be done using a Monte 
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Figure 2. Equation of state of quark nuclear matter. The solid line is for matter composed 
of nucleons only and the dashed line is for matter composed by nuclcons and quarks. 



Carlo integrator. For our purposes here, in this initial investigation we make 
some approximations. Initially we neglect the lower component of the Dirac 
spinor. This does not seem to be a too drastic approximation, since 7 in Eq. (4) 
is a small quantity. In this approximation, one obtains 

Z- 1 = l+ f d 3 hd 3 k 2 d 3 k 3 \%(h,k 2 , fc 3 )l 2 W kl ~ y 1 , 2 . , (26) 
J AE(p,ki,k 2 ,k 3 ) 

with F(q) given by 

F(q) = f d 3 kAV(k)e- k2 / x2 (e^ x2 - l), (27) 

where 3> p (ki, k 2 , k 3 ) is the three-quark wave function of the nucleon with cm. 
momentum p (see Ref. 4 ), AE(p, k\, k 2 , ks) is the difference between of the 
energies of the three unbound quarks and of the three quarks bound in the 
potential, and AV(k) is the Fourier transform of AV(r), where 

AV(r) =ar(e-^ r ' 2 - l) . (28) 

This clearly shows that once p, = 0, i.e. the potential is density independent, 
one regains the original QMC model. 

Now, Eqs. (26) and (27) still require a lot of numerical work. We simplify 
them further by making two additional approximations. The first one consists 
in neglecting the momentum dependence of the energy denominator and the 
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second one is to use an average value for q 2 in F(q 2 ). Both approximations 
taken together seem not to be a bad approximation, since the energy denom- 
inator under the integral is dominated by low momenta. Now the problem 
consists in a single one dimensional integral that can easily be performed with 
a Gauss integration. 

In Fig. (2) we present the results for the energy per baryon number, E/B 
as a function of the ratio pb/po- The solid line in this figure is the result 
for the QMC model of the previous section. The dashed line shows that the 
deconfining of quarks leads to a softening of the equation of state. It would 
be interesting to investigate the consequences of this softening for neutron-star 
phenomenology Soft equations of state seem to be required to explain recent 
observational data of compact stellar objects. 

4. Conclusions 

We have generalized the QMC model to include quark dcconfincmcnt in mat- 
ter. The model is based on an effective quark-hadron Hamiltonian obtained via 
a mapping procedure from a relativistic microscopic quark Hamiltonian with 
a density dependent quark-quark interaction. The equation of state of QNM 
was obtained using the effective quark-hadron Hamiltonian. It was found that 
beyond a critical density, when quarks start to deconfine, the equation of state 
predicted for QNM is softer than the usual QMC equation of state. Implica- 
tions of this equation of state for the phenomenology of compact stellar objects 
were pointed out. 
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